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ABSTRACT 


The separated flow environment induced by underexpanded rocket plumes during boost phase 
of rocket vehicles has been investigated. A simple semi-empirical model for predicting the 
extent of separation was developed. This model offers considerable computational economy 
as compared to other schemes reported in the literature, and has been shown to be in good 
agreement with limited flight data. The unsteady pressure field in plume-induced separated 
regions was investigated. It was found that fluctuations differed from those for a rigid flare 
only at low frequencies. The major difference between plume-induced separation and 
flare-induced separation was shown to be an increase in shock oscillation distance for the 
plume case. The prediction schemes were applied to PRR shuttle launch configuration. It 
was found that fluctuating pressures from plume-induced separation are not as severe as for 
other fluctuating environments at the critical flight condition of maximum dynamic pressure. 
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NOMENCLATURE 


Sound speed 

Local thickness of separated region 

Plume drag 

Frequency 

Frequence corresponding to macroscale of spectrum 

Separation region response frequency 

Plume response frequency 

Boundary layer shape factor 

Separation length 

Vehicle length 

Mach number 

Pressure 

Engine chamber pressure 
Plume nose radius 
Vehicle radius 
Reynolds number 
Engine thrust 
Speed 

Convection velocity 
Free stream speed 
Streamwise coordinate 
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Ratio of specific heats 
Boundary layer thickness 

Boundary layer thickness immediately upstream 

Boundary layer displacement thickness in transformed coordinates 

Boundary layer momentum thickness in transformed coordinates 

Plume angle at nozzle lip 

Reattachment turning angle 

Separation angle 

Spectral density 

Condition at separation point 

Condition in attached turbulent boundary layer 

Condition in separated region 

Condition at shock wave location 

Average 


1.0 INTRODUCTION 


During atmospheric flight, the external surfaces of the space shuttle vehicles will be subjected 
to high intensity acoustic and fluctuating pressure environments. Wyle Laboratories is engaged 
in a research program under NASA-MSFC Contract NAS8-26919 to develop reliable methods 
of predicting critical acoustic and fluctuating pressure environment trends during space shuttle 
flight. In the course of this program, fluctuations associated with various attached and 
separated flow environments have been defined and quantitative prediction formulae developed 
(References 1-7). Among the most severe fluctuating environments are those associated with 
compression induced separated flow and shock oscillation. These are associated with protu- 
berances and flares in the vehicle geometry, and, as described below, exhaust plumes. 

During the boost phase of atmospheric flight, engine exhaust material forms an effective 
afterbody. When the exhaust is highly underexpanded (nozzle exit pressure much greater than 
free-stream pressure) the resultant plume can be significantly larger than the vehicle itself. 
Under these conditions, compression induced separation may exist over a significant area of 
the vehicle due to plume interference. When the Reynolds number is high enough so that the 
flow field may be treated as inviscid flow plus boundary layer (with strong interactions occurring 
only locally at separation points), plume effects may be treated by considering the plume to be 
an afterbody whose surface corresponds to the dividing streamline between exhaust material 
and the ambient. Calculation of the flow field is more difficult than for a rigid body because 
the plume shape is implicitly dependent on the pressure field. Because the equivalent plume 
"body" is a dividing streamline and may respond to pressure fluctuations, the fluctuating flow 
environment cannot be assumed a priori to correspond exactly to that of a rigid afterbody. 

In this report, the fluctuating pressures associated with plume induced separation are examined. 
The basic properties of supersonic flow separation and of highly underexpanded rocket plumes 
are first discussed, and an approximate model for separation length and steady properties is 
presented. Although more elaborate models may provide more accurate predictions of 
separation lengths, the simple model employed clearly exhibits the overall scaling properties 
of plume-induced separation. The flow properties deduced from this model are then applied to 
the known properties of fluctuating pressures associated with separation by rigid flared bodies. 
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Anticipated departures from these properties due to the non-rigid nature of the plume are 
identified. These include fluctuations of the plume itself due to combustion instabilities. 
Finally, predictions of boost phase plume-induced fluctuating pressures are made for a 
preliminary shuttle configuration and mission. 
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2.0 PLUME-INDUCED SEPARATED FLOW FIELDS 


The basic flow configuration of interest ?s illustrated in Figure 1 . There are three distinct 
regions: exhaust flow, outer separated flow, and recirculation region. The turning angles of 
the outer flow at the separation point and reattachment point are governed by viscous interaction: 
the turning angle must be consistent with the increase in boundary layer displacement thickness. 
The locations of the slip lines are determined by pressure matching. The reattachment point is 
located such that the pressure rise across the inner plume shock equals the pressure rise across 
the reattachment shock. Boger, Rosenbaum and Reeves (Reference 8) and Fong (Reference 9) 
have developed computer programs to calculate the separation induced flow field for an 
initially laminar boundary layer. The individual components of the flow field shown in 
Figure I are discussed in detail in the following subsections. 

2.1 Viscous Separating and Reattaching Flows 

The problem of separation of supersonic flows is an important one for high speed vehicles, and 
has been the subject of a number of theoretical and experimental investigations. Figures 2 and 
3 show two supersonic configurations leading to separated flow: an impinging shock wave in 
Figure 2, and a compression corner in Figure 3. An important feature of both flows is that 
information of compression is propagated upstream through the boundary layer. Without a 
boundary layer, the separation region would degenerate into a point — the shock impingement 
point in Figure 2 or the corner in Figure 3. It is clear, therefore, that compression-induced 
separation is a viscous interaction phenomenon. 

Various theoretical treatments have shown that boundary layers may be supercritical or 
subcritical . In a subcritical boundary layer, disturbances are communicated upstream; in a 
supercritical layer they are not. There is a direct analogy to subsonic and supersonic flow. 

In the case of a supercritical boundary layer, communication of a pressure rise can be achieved 
only through a "jump" to subcritical, just as with a shock wave. After reattachment, the 
layer goes through a neck region where a transition from subcritical to supercritical takes 
place. Crocco and Lees (Reference 10) showed that the correct solution for reattachment is 
one which passes smoothly through this critical point. 
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Theoretical approaches have generally been to use the integral form of the boundary layer 
equations. Early studies were for simple methods, such as Karman-Pohlhausen, with more 
sophisticated methods being gradually adopted. Theoretical development was also first for 
laminar boundary layers (Reference 11), being modified later to include turbulent layers. The 
most recent theoretical model of turbulent separation is due to Hunter and Reeves (Reference 12) 
and Todisco and Reeves (Reference 13). In this model, the laminar sublayer and skin friction 
are neglected, so that the .low is wake-like. This appears to be a reasonable approximation, 
as skin friction does vanish at separation. Boundary layer behavior at distances upstream of 
separation does not affect this model as virtually all attached supersonic turbulent boundary 
layers are supercritical (Reference 14). 

The model developed by Todisco and Reeves and extended by Hunter and Reeves applies to 

the separated flow configuration of Figure 3. A set of supercritical-subcritical jump conditions 

(analogous to normal shock relations) were found for the separation point. After separation, 

development of the free shear layer is calculated by numerical integration of the moment 

integra. form of ihe boundary layer equations. Shear layer profiles are expressed as a one 

parameter function of chape factor H = 0. /8* , where 0. and 5* are momentum and 

11 11 

displacement thickness, respectively, in transformed coordinates. After reattachment, the 
pressure is assumed to be the inviscid wedge pressure corresponding to the wedge angle. The 
wedge angle is adjusted until the solution passes smoothly through the critical point. 

Figure 4 shows a typical result obtained by Todisco and Reeves. Flow conditions are noted on 
the figure. A separation shock with pressure ratio of about 2 (this does not increase markedly, 
even at hypersonic Mach numbers) is followed by a rise to plateau pressure. At the corner, 
pressure rises more sharply, approaching ramp pressure downstream of the critical point. 

The variation of plateau pressure with Mach number is shown in Figure 5, along with some 
experimental data. According to the Todisco and Reeves theory, plateau pressure is independent 
of ramp angle. This has been questioned - at least for the hypersonic case - by Elfstrom 
(Reference 15). Figure 6(a) shows plateau pressures measured by Elfstrom added to Figure 5. 
Figure 6(b) shows this Jata normalized with respect to inviscid ramp pressure, with apparently 
excellent collapse. However, Reeves (Reference 16) has pointed out that all of Elfstrom 's 
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measurements are far separation lengths no greater than 4 boundary layer thicknesses; at 
such small lengths, pressure at the corner has not reached plateau conditions. Reeves presents 
plots (Figure 7) which show reasonable agreement when p/p^ ; s collapsed according to 
Ax/S . Most available data appears to be for separation lengths l / 6 < 8 , at which point 
plateau conditions ore still not reached. Reeves pointed out that experiments capable of 
providing data for i^/6 at least 15 should be performed, and would settle this question. 

For the present, it appears that Reeves' interpretation of the data is supported by the Todisco- 
Reeves theory, and will therefore be adopted in this report. 

An interesting feature of turbulent separation is that Reynolds number dependence is very weak, 
becoming negligible for Reynolds number greater than two or three times the value at transition 
(Reference 17). This would appear at first to be peculiar for a viscous interaction phenomenon. 
However, it should be noted that at high enough Reynolds numbers the inertial subrange will 
dominate the turbulence vectrum, so that viscosity itself has a less important role. This is 
consistent with, for examp. 3, the empirical correlation for turbulent boundary layer thickness 
(Reference 18) 



( 1 ) 


which becomes 5 = 0.01 x at large R g . This is a considerable simplification over laminar 
separation, in which Reynolds number remains an important parameter. Separation of a 
turbulent supersonic boundary layer should depend only on local boundary layer properties and 
Mach number. For zero pressure gradient, only Mach number remains as an important 
parameter. For the purpose of this report, therefore, separation plateau pressure will be taken 
as that given by the Todisco-Reeves theory, shown in Figure 5, and separation angle based 
on tangent wedge for this pressure. 


5 




* 


d 


2.2 Plume Flow Field 


Figure 8 illustrates the inviscid flow field of a supersonic highly underexpanded rocket plume. 
Region I is free-stream flow. Flow in Region IV is exhaust flow, and is the same flow as 
would exist in a vacuum. Flow in Regions I and IV must both turn in order to match pressure 
and direction at the contact surface; since both flows are supersonic, they will experience 
shock waves. Flow within Region II has been turned outward by the outer plume shock, with 
the flow in Regions I and II being the same as for a body with shape corresponding to the 
interface. Flow in Region IV is turned by the barrel shock so that flow in Region III is more 
in the axial direction. Since pressure in Region IV diminishes rapidly downstream as the 
vacuum flow expands radially, the barrel shock must increase in strength with downstream 
distance. At the nozzle lip, where engine exhaust simply expands isentropically until it 
matches pressure at the interface, the barrel shock has zero strength; far downstream it 
becomes stronger, approaching the hypersonic limit. At some point downstream the barrel 
shock terminates at a normal shock, generally referred to as the Mach disc. A series of "cells" 
similar to the Mach diamonds of a slightly overexpanded jet, may exist beyond this; however, 
for the present problem only the forward portion of this first cell Is of interest. 

In the case of a plume in a quiescent ambient. Region II does not exist and Region I has zero 
velocity and constant pressure. Calculation of the plume flow field is then a straightforward 
matter using the method of characteristics, and has been do <e by a number of investiaotors, 
including References 19 and 20. When the ambient is flowing, the pressure becomes 
a function of the interface shape. If conditions are such that the outer shock is attached to 
the nozzle lip and the flow Is everywhere supersonic, a method of characteristics calculation 
should present no difficulties, although the calculation would be more complex because 
exhaust material and ambient have different properties. 

At very high altitudes, the initial exhaust expansion at the nozzle lip is great enough so 
that an attached outer shock cannot exist (Reference 21 ); in this case the bow shock is 
detached and Region II corresponds to the flow about a blunt body. Supersonic blunt body 
flow calculations are difficult enough with a known body shape. The combination of a state 
of the art blunt body computer program with a method of characteristics program presents 
formidable difficulties and has not yet been done. 


Because a computation of high speed plume flow fields using the most advanced techniques 
presents great difficulties, calculations to date have used various approximate methods. 
Boynton (Reference 22) has employed a Lagrangian (stream tube) finite difference method. 

It is not ideal for supersonic flows, as Boynton points out, because it does not implicitly 
contain the wave nature of such flows and shock waves present special problems. Its main 
advantage is providing some answer where no other exists. Another approach is to use 
Newtonian flow to locate the dividing streamline (References 23 and 24). The primary 
disadvantage of this method is that the inner shock is very weak near the nose of the plume, 
so that Newtonian flow is not a very good approximation for the nose portion of Region HI. 


Although detailed calculations of plume flow fields have not yer been entirely successful, 
certain scaling properties have been well established. Jarvinen and Hill (Reference 25) 
have developed a semi-empirical "universal plume profile" which incorporates these scaling 
laws. Figure 9 shows the Jarvinen-Hill model. Dimensions in the axial direction scale as 
^ ^ ' wliere T = rockef ^rust and i: freestream dynamic pressure. Dimensions 
In the lateral direction scale as ((TD)5/q^)s , where D= plume drag. Plume drag is definec 
as the vacuum thrust which would result if the exhaust were fully expanded, minus the actual 
vacuum thrust. The physical interpretation of this quantity may be seen if the interface is 
imagined to be a nozzle wall; D is simply the total aerodynamic force exerted on the plume 
material. The universal plume profile shown in Figure 9 is based on experimental data, 

method of characteristics computation and Newtonian flow computations, with the nose region 
based on a blasf wave analogy. 


The Jarvinen-Hill plume profile provides a good working model except for the very nose, 
where detail necessarily neglected by such a model Is important. The parobolic nose 
resulting from the blast wave analogy implies a 90° turning angle of exhaust at the nozzle 
lip; this is generally not correct. The initial angle of the plume is, however, simple to 
compute. Reis, et al. (Reference 20) noted that, for plumes into a quiescent ambient, only 
one family of characters! ics entered the computation very near the nozzle lip. This means 
that the inner flow for a very short distance from the lip can be determined by a Prandtl-Meyer 
expansion. If the outer flow pressure Is obtained by Newtonian flow or tangent cone, the 
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initial angle of the contact surface is easily obtained. The universal plume profile can then 
be modified to include the actual angle for each set of flight conditions. 

2*3 Plume-Induced Separation 

Boger, Rosenbaum and Reeves (Reference 8 ) (hereinafter called BRR) and Fong (Reference 9 ) 
have recently published schemes for computing plume-induced separation. The BRR 
scheme is described in Section 2.3.1 ; Fong's scheme is roughly similar but more 
elaborate, including the shear layer at the plume interface as well us the separation line. A 
simplified scheme is presented in Section 2.3.2. 

2.3.1 Sep aration Scheme of Boger, Rosenbaum and Reeves — The BRR scheme follows 
the basic elements outlined in Section 2.0: separation and reattachment are governed by viscous 
interaction, with the plume boundary determined by an inviscid calculation. The numerical 
scheme follows four sfeps: 

1 ) A separation point is assumed. Separation angle, plateau pressure and 
Mach number corresponding to plateau pressure are determined. The 
example used in Reference 8 was for an initially laminar attached 
boundary layer; plateau pressure was computed using a relation from 
Reference 26; angle and Mach number were obtained from this by 
tangent-cone. 

2) Using the plateau pressure, the plume boundary is computed by the 
method of characteristics. Because plateau pressure is constant, this 
was done by a scheme developed for computing plumes in a quiescent 
ambient. Neglecting the viscous layer at the plume boundary in deter- 
mining the boundary is justified by a series of wind tunnel experiments 
reported in Reference 8 . The plume boundary, separation angle, and 
assumed separation point give a preliminary reattachment point on the 
plume boundary. 
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3) The turning angle at reattachment is determined from the viscous interaction 
theory of Todisco and Reeves (Reference 13). For a given plateau Mach 
number and shape factor H = 0./8* , the turning angle is uniquely deter- 
mined. Figure 10 shows the conditions used by BRR. At the separation point 
H = H s , so that H/H s = 1 ; as the shear layer develops, H/H decreases, 
approaching zero for a fully developed free shear layer. 

4) The pressure after reattachment is known from the turning angle determined 
in Step 3. The point on the plume boundary which gives the same pressure 
after turning to the same final angle is located. If this point does not agree 
with the point from Steps 1 and 2, a new separation point consistent with 
this is chosen. The four steps are iterated until the solution converges. 

The BRR scheme as described here requires recomputation of the plume boundary at each iteration 
only if plateau conditions are dependent on separation point. For a smooth body with zero 
pressure gradient at the separation point, this means a dependence on Reynolds number. The 
laminar example presented in Reference 8 has a Reynolds number dependence. For a turbulent 
attached boundary layer at very large Reynolds number, as described in Section 2.1, there is 
no Reynolds number dependence. Iteration of the separation point is still required, however, 
because the reattachment conditions depend on the length of the shear layer. 

— * . 2 A PP ro * imq _te Model for Plume-Induced Separation - The BRR scheme appears 
to be a reasonable method for computing plume-induced separation. However, there are 
certain approximations in calculating the separated flow. The most significant is assuming 
constant pressure in the plateau region. As seen in Figure 4, this is clearly not the case: 
pressure varies somewhat after separation, and considerably after the corner. This latter is 
crucial, as it is the assumption of constant pressure that enables a characteristics solution 
of the plume flow. A proper calculation, taking the pressure gradient into account would be 
of the same difficulty as the hypersonic plume with detached bow shock discussed earlier. 

Other approximations, such as assuming a constant separation angle, are not as significant 
as this. 
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Since the boundary conditions for the method of characteristics plume calculation are not 
correct, it is worth seeing if further approximations are possible so as to eliminate entirely 
the need for it. As this is the most complex part of the BRR scheme, it would greatly 
simplify computation. Toward this end, the following assumptions are made: 

1 ) As with BRR, pressure within the separation region is assumed constant. 

2) The plume boundary adjacent to the separation region is adequately 
described by a Prandtl -Meyer expansion to match plateau pressure. 

The contact surface is thus conical in this region. This is the same 
model as used in Reference 21 to predict detachment of the outer 
plume shock. 

3) Outside of the separated region, the Jarvinen-Hill universal plume 

profile is adequate to describe the contact surface. In the nose 

region, the surface is a sphere with radius r . 

P 

With these assumptions in hand, the plume-induced separation flow model is shown in Figure 11 
The separation angle 0 and plateau pressure and Mach number are given by the Todisco— 
Reeves theory, 9 p! is found by expanding exhaust flow to match plateau pressure, and 
r p = 0.2 (D/q^ )? is from the Jarvinen-Hill plume model. The reattachment turning angle, 
the separation point and center of spherical interface are such that the reattachment angle 

(0 sep + 0 reat^ ' S tan 9 enf to the s P here at * he reattachment point. The separation length is 
then given by 

* = [ V COi <( W + 9 r.a.> - 'v ] ( C °' e !ep - “* 9 p| ) < 

For typical vehicles of interest here, the boundary layer is fully developed turbulent. 
Separation pressure is therefore given by Figure 5, with Mach number and angle determined 
from the wedge flow charts of Reference 27. Separation values at M < 2 were extrapolated. 
The values of M and 0 used are shown in Fiaure 1 2. 

sen cpn ^ * 



Because the vehicles of interest, and consequently separation lengths, are very large, it is 
expected that the separated shear layer at reattachment will be close to fully developed. 

Therefore, reattachment conditions are taken as given by the curve for H/H =0.11 on 

Figure 10. This curve is chosen because it is the most developed case presented in Reference 
8 . Figure 13 shows this curve extrapolated down to Mach numbers of interest here. 

A computer program, described in Appendix A, was written to compute separation conditions 
versus time for any given flight trajectory. The program includes fluctuating conditions, develop- 
ed in Section 3.0, as well as steady separation conditions discussed thus far. 

As an example of this simplified calculation, separation conditions were calculated for the 
Saturn V S-l , using a nominal Apollo launch trajectory (Reference 28). Trajectory parameters 
are shown in Figure 14. Figure 15 shows the predicted separation point versus time after 
lift off. Data (Reference 29) obtained from a microphone at station 757 on an 
Apollo flight indicated large fluctuating pressure of the type associated with shock wave 
oscillation from approximately 118 seconds to approximately 130 seconds. This range is 
indicated on Figure 15; the separation point passes station 757 somewhere during this time 
period. Since the shock wave precedes the separation point by a small amount, the separation 
point most likely passes this station somewhat after the middle of this time range. Agreement 
betv/een this data point and the prediction is good. The prediction method presented in this 
section for plume-induced separation is therefore quite adequate to use as a framework for 
estimating fluctuating pressure environments. 
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3.0 FLUCTUATIONS IN PLUME-INDUCED SEPARATED ENVIRONMENTS 

3.1 Summary of Fluctuating Pressure Environment in a Compression Corner 

Figure 16, from Reference 30, shows the static and fluctuating pressure environments ahead of 
a compression corner. The fluctuating pressure may be divided into three parts: incoming 
boundary layer, shock wave oscillation, and homogeneous separated flow. The spectra 
associated with each of these are shown in Figure 17. Reference 4 reviews available data and 
provides empirical prediction schemes for these environments. The amplitudes and spectra 
shown in Figures 16 and 17 are representative of ‘wo dimensional and axisymmetric compression 
induced separation; the scaling parameters used in these figures were found in Reference 4 . 
to give best collapse of experimental data. Application of prediction schemes to specific 
vehicle configurations is straightforward, and shows excellent agreement with wind tunnel 
measurements (Reference 7 ). 

Although empirical predictions of fluctuating properties can be made with some confidence, 
the fundamental mechanism involved are not fully understood for all parts of the flow. The 
following properties, relative to assessing the effect of the free plume boundary, are established: 

3.1 .1 Turbulent Boundary Layer - This is the best understood part of the flow. Since 
supersonic turbulent boundary layers are generally supercritical, boundary layer properties 
upstream of the separation point will not be affected by the plume boundary. 

3.1.2 Sh ock Wave Oscillation — An analytical model has been established linking 
the motion of the shock wave to fluctuations in the incoming turbulent boundary layer 
(Reference 6 ). The shock wave is convected by velocity fluctuations in the boundary layer. 
Since boundary layer fluctuations may be treated as a stochastic process, the shock wave 
executes a random walk about its mean location. Mean flow properties seek to return the 
shock to its mean location, so that displacement is limited. At large times, the shock motion 
is a stationary random process with mean square displacement dependent on both boundary layer 
fluctuations and the mean flow restoring mechanism. The macroscale of the displacement is 
governed by the mean flow restoring mechanism. Conditions within the separation region may 
therefore influence the shock motion. For example, slowly varying fluctuating pressure in the 
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exhaust flow would cause the separation length, hence shock location, to fluctuate. (For 
exhaust oscillations slow compared to response time of the separated flow, shock motion may 
be calculated as quasi -steady.) For a steady exhaust flow, it is expected that the nature of 
the shock motion (driven by the boundary layer) will not change, but that the total 
excursion may. 

3.1 .3 Separated Flow Region — No quantitative model has yet been developed for 
this region. Even mean flow conditions are not yet fully established (Reference 15). However, 
the fluctuating flow is qualitatively understood and arguments may be put forth as to the 
mechanisms of the fluctuations. 

The mean flow in a compression corner has been partly described in Section 2.1 , and is shown 
In Figure 3. A key feature of the flow for interpreting fluctuations is the region beneath 
the separation streamline. In the discussion of Section 2.1 , this was treated essentially as a 
"dead water" region. As with mosr "dead water" regions, it is actually a recirculation zone, 
with trapped eddies. For two-dimensional and axisymmetric flows, there is a single large 
eddy in the separation region. This eddy is a key factor in the fluctuating pressure field, 
expecially at low frequencies. With a single trapped eddy, two-dimensional and axisymmetric 
flows should have similar fluctuating pressures; this is borne out by Figures 16 and 17. 

Separated flow induced by three-dimensional protuberances, on the other hand, exhibits 
different (and larger) fluctuating pressure levels (Figure 18). This is clearly a consequence 
of the multiple trapped eddy structure, identified by Robertson's oil flow experiments 
(Reference 2). For the present problem, only twc-dimensional and axisymmetric flows are of 
interest. However, it should be noted that the turbulent boundary layer and shock oscillation 
fluctuations are essentially the same for protuberances as for two-dimensional flow. This is in 
agreement with the conclusions of Sections 3.1 .1 and 3.1 .2. 

Interpretation of fluctuations In the separation region requires knowledge of narrow band 
convection velocities. Figure 19 shows convection velocities presented by Coe and Rechtlen 
(Reference 30). Convection velocities presented by Chyu and Hanley (Reference 31 ) and 
Rechtien (Reference 32) are essentially the same. 
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Reference 32 does indicate negative convection velocities at very low frequencies (approxi- 
mately f 6/ Uao <0.03), but there is considerable scatter in that range. Because of the scatter, 
and because low frequency random data is inherently less accurate than higher frequencies, it 
is not clear whether or not negative convection velocities really exist. The phase angle of the 
cross -correlation function is also near zero in this range, which makes accurate calculation 
of convection speed difficult. 

The convection velocities shown in Figure 19 vary from a minimum of u /u ss 0.2 at 

fS / U oo ~ 0 * 06 Pl^au ^/u^ « 0.8 at fS/u^O.8. Shown for comparison with the 

plateau are convection velocities obtained by Bull (Reference 33) for attached turbulent 
boundary layers. Agreement with the plateau above fS/u^-0. 8 is good. Noting that eddies 
of size 8 convected at U £ = 0.8 u^ correspond to fS/u^O.8, it may be concluded that 
fluctuations at high frequencies, fS/u^O.8 , are associated with the free shear layer. 

Rechtein (Reference 32) provides qualitative arguments that pressure fluctuations in this 
frequency range are transmitted to the wall by the "eddy Mach wave" phenomenon (Reference 
34). For the present ourposes, however, it is sufficient to only identify this frequency range 
with the free shear layer. It may then be concluded that this part of the fluctuations will not 
be affected by having a non-r;gid plume boundary downstream. 

For f8/ Ua <0.8, attached boundary layer convection velocities slowly rise, approaching 

U qo ° S f5/u aT° ( Referenc e 33), while convection speed in the separation corner decreases. 
This difference is because low frequency boundary layer noise is identified with free stream 
air entrained at the edge of the layer, while the free shear layer et., rains nearly still air from 
the separation bubble under it. 


Convection velocities for f8 /u oo < 0.8 can be interpreted in terms of a simple entrainment 
model. The free shear layer is initially composed of eddies of size 8, convected at u = 0.8 u 
Essentially still air beneath it is entrained, with the vorticity of the original layer distributed °° 
over the entire layer. If the thickness of the layer is now d, eddy sizes are also of size d. 

The convection speed of the larger eddies is 


u 

c 


0.8 u i- 
od d 


(3) 
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This represents conservation of momentum of the original layer. The frequency associated with 

a series of eddies of size d moving at speed u is 

c 

u 


so that 


f6 = _6_ 

u d u 

00 00 



(4) 


Figure 20 shov/s a model convection speed, based on boundary layer convection for f5/u £0.8 
and Equations (3) and (4) for fS/u^0.8. Also shown are convection speeds from Figure 19. 
Agreement of this simple entrainment model is quite good for f6/u larger than about 0.1 
which corresponds to d/6 * 3. In the experiments of Reference 30, step height was not more 
than a few boundary layer thicknesses. Since the entrained layer cannot be thicker than 
approximately the step height (and must be somewhat less because there must be a reveree flow 
near the wall), it is not expected that this entrainment model would apply at lower frequencies 
than this. 

The agreement between convection speeds for the entrainment model and the data indicates 

that fluctuations in this frequency range are also generated by the shear layer. The plume 

boundary should therefore have no major effect on fluctuations above f6/u *01 

oo 

The convection velocity has a minimum of about 0.2 u^ at fb/u^ = 0.06, then rises sharply 
below that. It Is in the region of this rise that Rechtein (Reference 32) obtained widely 
scattered positive and negative convection velocities. Since a solenoidal pulsation would 
result in infinite positive and negative computed convection speeds, it may be concluded 
that the region fS/u^ 0.06 is dominated by large scale fluctuations of the whole separation 
bubble. This conclusion may be supported by considering the acoustic response frequency of 
the separation bubble. If the separation length is 106 (experimental measurements cited here 
are generally less) and M = 2, then the fundamental response of a cavity of equal size 
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corresponds to fS/u^- 0.05. This is approximately where large scale fluctuations begin to 
dominate over the convected eddies. In this range, it is expected that there will be coupling 
between the plume boundary and fluctuating separation pressure. 

3.2 Coupling Between Plume Boundary and Separation Fluctuations 

The discussion of Section 3.1 leads to the conclusion that above a certain frequency separated 
flow fluctuations are dependent primarily on the incoming turbulence. Below this frequency, 
the geometry of the mean flow (i.e., the actual dimensions of the separation region) has a 
direct effect on fluctuations. It is in this frequency range that flow over a plume may differ 
from flow over a rigid flare. 


Because the feedback mechanism through the plume interface involves motion of the plume, 
it is necessary to assess the response of the plume to pressure fluctuations. There will be some 
frequency (based on characteristic flow time of the plume) below which the plume will exhibit 
direct response, following pressure fluctuations quasi-statically. Coupling of motion in this 
regime will be relatively straightforward. At higher frequencies, response of the plume will 
lag somewhat, with direct feedback being less as fluctuation frequency increases. 

In Section 3.2.1 below, the various characteristic frequencies are summarized and the 
frequency range of separated flow fluctuation/plume boundary coupling identified. The 
quantitative effect on fluctuating pressure environments are then estimated In Section 3.2.2. 

3.2.1 Frequency Range of Fluctuation/Plume Coupling - The "fundamental" 

frequency of a separated region of length l is 

s 

a 

r 00 

sep 2 1 (5) 

s 

where it is assumed that sound speed within the separated region is a . This is the first 

acoustic mode of a cavity of size I . At frequencies less than f ® pressure fluctuations 
are in phase. This corresponds to the infinite convection speeds discussed in Section 3.1 .3. 
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Casting f in the form f 8 / u , 
sep 0 go 


a 5 

f 8 /u = -2° -± 

sep o oo 21 u 

s ao 



(6) 


The value of f given by Equation (6) corresponds to the lower frequency limit 

of convected fluctuations, as discussed in Section 3.1 .3. It was concluded in that section 
that at frequencies greater than f the presence of a non-rigid plume boundary does not 
influence separated flow fluctuations. 

For the plume interface to follow pressure fluctuations quasi -statically, the plume flow time 
must be greater than the fluctuation flow time. The plume flow time is given by the length 
of plume interface within the separated region, divided by flow speed. For the present 
purposes, the length may be taken as l tan 0 and the speed as u . For typical 
separation angles tan 0 ss 1/4 , so that the response frequency of the plume is 

4u 

f = — i 
pl l 

s 


Casting in the form f6 /u . 

0 CD 


f ,6/u 
pl 0 00 



00 


8 

0 


(7) 
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Comparing Equations (6) and (7), 


u i 

f . = 8 -1 JL 

pi u M 

00 00 

so that f , is generally an order of magnitude larger than f . This me^ns that at 
p sep 

frequencies where pressure fluctuations are coupled to the plume motion, the plume boundary 

response is quasi-steady. This is an important finding, and much simplifies calculation of the 
response. 

The frequencies given by Equations (6) and (7) must be compared to characteristic frequencies 
of the fluctuating flow. The most important frequencies are those corresponding to the integral 
scales of the fluctuation spectra. The bulk of the overall fluctuating pressure level is 
identified with a frequency band about this value. 

The integral time scale of a fluctuating environment is defined in Reference 6; the corres- 
ponding frequency is the reciprical of the time. For the spectra shown in Figure 17, the 
characteristic frequencies for the boundary layer, separated region, and shock wave 
oscillation are 


sep 


( 8 ) 
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6 /u 
0 00 
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sep 
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00 


f 6 / u 
o o oo 
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0.53 


0.145 


0.04 


(9a) 

(9b) 

(9c) 
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Figure 21 shows Equations (6), (7) and (9) as a function of t /S . Also indicated is £ /L , 

so s 

separation length divided by total vehicle length. This was obtained by assuming 6 o is given 

by the large Reynolds number limit of Equation (1 ), and noting that distance from the front 

of the vehicle is L - l .so that 

s 


S 0 = 10 " 2 ( L -'V 00) 

Comparison of the various curves in Figure 21 leads to the following conclusions: 

• Except for small separation lengths (£ /6 < 1 at M = 2, for example), 

SO CD r ' 

overall fluctuating pressure level in the separation region is not affected, 

• For ^/6 £ 10, shock wave motion may be strongly coupled to plume 

motion. It should be noted that shock oscillation is not a homogeneous 
environment, so that important low frequency effects may occur at larger 
separation lengths as well. However, although these effects may be 
important, it is not expected that overall levels will change significantly. 

• Boundary layer fluctuations upstream of separation, and high frequency 
fluctuations in the separation region, will not be affected. 

3.2.2 Amplitude of Coupling - For all important frequencies, plume response to 
separation fluctuations will be quasi-static . The plume size changes in response to separation 

pressure fluctuations. The separation point and shock wave then move in accordance with the 
instantaneous plume size. 

The mean square pressure fluctuation in the separation region Is 


< P ? > 



(f) df 


(in 
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Feedback to the separated region is only for f 6 / u < 1/M 5 // . f he pressure 

u 00 00 0 s r 

fluctuation in this frequency range is 


A < p 2 > 



(12) 


Equation (12) has been numerically integrated, using the analytic representation of 0 

1 1 sep 

presented in Reference 4. The quantity A<p 2 > 5 /< p 2 >2 is p l ofted as a function of 

M ^ seD / 8 0 in Figure 22. Figure 23 shows A<p 2 >Vq as a function of / /6 for 

“ u 00 sep 0 

several Mach numbers, based on Figure 22 and the Mach number dependence of < p 2 > * /q 

00 

as presented in Reference 4. The quantity A < p 2 > Vq may be considered to be the 

00 

feedback fluctuating pressure coefficient", since it represents that part of the fluctuating 
pressure which can couple directly to the plume. 

The response of the plume to the feedback fluctuating pressure is straightforward to compute. 
Referring to Figure 1 1 , 0 ^ , 0 ^ and r p are functions of the external flow only, so that 
the only change in geometry is through changes in 0 . If 0 changes by an amount A0 

Pi PI ' n| ' 

the separation length changes by 


A ^S wTT 


sep 


cot 0 


tan A 0 


pl 


Pi 


1 + cot 2 0 


pl 


1 + tan A0 . cot 0 . 

pl pl 


( 13 ) 


For small changes in A0 p( , such that tan A0 p[ « 1 , Equation (13) may be represented as 


A f 
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6.1 
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(14) 
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where 
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(15) 


Root mean square change in separation length due to change in plateau pressure is given by 



3 . d9 p' ( A<p2> ) j 

d9 pl d ( p/ ' , CD ) 


(16) 


where d0^ / d (p/q^) must be obtained for the plume flow field. Since 0^ ( is obtained 
by a Prandtl-Meyer expansion of the exhaust flow to match plateau pressure, it is straight- 
forward to obtain 


d 0 


d p/q 
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where M is the Mach number of the exhaust flow after expanding to 0 

pl 

obtained by differentiating Equations (2.31 ) and (4.21b) of Reference 36. 


Equation (1 7) is 


All parameters needed for Equations (15) and (17) are available from the calculation of 
separation length. Figure 24 shows d/ 8 /d(p/q ) versus Time After Lift-Off for the 

Saturn V, nominal Apollo trajectory. Using the feedback fluctuating pressure coefficient 
shown in igure 22, the change in separation length has been calculated. Figure 25 shows 
the increase in root mean square separation length due to the flexible plume boundary; 
shown for comparison are boundary layer thickness and root mean square shock oscillation 

displacement for an equivalent rigid flare. The net rms shock motion will be no greater than 

1 1 

< x . w > * < Af ‘ > (this represents perfect coherence between x and A/). Fiaure 26 

3 sw 3 


t 
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o 1 

shows < A i 2 )' /< x 2 > 2 for the Saturn V. The Increase in rms shock wave displace- 
ment is less than 23 percent. Fluctuating pressure at the shock wave is proportional 
to displacement (Reference 6); the rms fluctuating pressure level at the shock location is 
therefore increased from 0.068 q (Reference 6) to 0.0836 q . On the decibel scale, this 

CD CD 

is an increase of 2 dB. 

Other than the slightly increased fluctuating level at the shock location, the change in 
separation length is not enough to significantly alter the fluctuating pressure patterns. It 
may therefore be concluded that, for the example presented, there is little difference 
between fluctuations associated with flare-induced separation and plume-induced separation. 
The only significant difference is a slight increase in shock oscillation distance and fluctuating 
level. 

The analysis so far has assumed that the exhaust flow itself is steady. The question of 
unsteady exhaust is treated in the next section. What may be concluded here Is that the 
fluctuating pressure due to plume-induced separation does not have an important coupling 
effect itself with the plume boundary. Homogeneous fluctuating pressure environments may 
therefore be superposed on fluctuating pressures associated with unsteady plume effects 
discussed in the next section. Pressure fluctuations due to the oscillating shock wave, which 
is an inhomogeneous environment, does not simply superpose; however, once the displacement 
is known the shock oscillation pressure field can be determined from the theory of Reference 6. 

3,3 Plume Driven Fluctuations 

The discussion thus far has assumed that the engine exhaust flow is steady, with turbulent 
fluctuations present only in the external separated flow. For real rocket engines this is not 
the case. At the high pressures and Reynolds numbers involved, flow within the combustion 
chamber is turbulent. Further, the combustion process, with vaporizing fuel droplets burning, 
is not uniform. Although the nature of fluctuations within a rocket engine can be of several 
types (Reference 35), for the present purposes only the effect of a solenoidal pulsing of 
combustion chamber pressure will be evaluated. 


22 


3.3.1 Low Frequency Fluctuations: Quasi-Static Response — If all other parameters 
are held fixed, thrust and exit pressure of a rocket engine are proportional to chamber 
pressure. For chamber pressure fluctuations of frequency less than f as shown in Figure 21 , 
the fluctuating plume shape may be computed quasi-staticaiiy, with T/< T > = p /< p > = 
P C /^P C ^ • Referring to Figure 11, the model plume shape will experience changes in 9 
and r p due to thrust fluctuations. The angles 0 and 0^ are unchanged. The plateau 
pressure is assumed to be unchanged; this assumption is correct for f < f , when separation 
region response is also quasi-static . 

For small changes in 6^ and r^ , the change in separation length is: 
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where is the same as dt^/6 6 . given in Equation (15), and 
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For small fluctuations in chamber pressure and thrust, A 0 and Ar are 
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Note that any appearance of AT/< T > in Equation (20) can be replaced by Ap /< p > . 

c c 

The plume radius used in the present model, described in Section 2.3.2, is 


r 

P 



( 21 ) 


Plume drag D is proportional to thrust T; thus, for small fluctuations in T, 

3r i 

P_ = I r 

3T/< T > 2 P 


( 22 ) 


The change in plume angle 0^ is governed by matching the Prandtl-Meyer expansion of the 

fluctuating exhaust flow to the constant plateau pressure. Derivation of d0 /dp // D ^ 

pl r c /N ' 

is similar to derivation of Equation (17), except that Equation (2.31 ) of Reference 38 is 
differential with respect to p^ rather than p . The result is 
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where M is the Mach number of the exhaust flow after expanding to 0 

P 1 * 

Combining liquations (18) through (23), 
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where 
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where 9^/90^ and ai $ /9r are given by Equations (15) and (19). 

Figure 28 shows Equation (25) for the Saturn V. Little data is available on fluctuating thrust 
of the F-l engine; howevet, if thrust fluctuates by 1%, separation length would vary by 
approximately 0.5 meter. As with the fluctuations discussed in Section 3.2, variations of 
separation length of this order do not significantly affect the homogeneous separated flow 
fluctuations. Fluctuations due to the increased motion of the shock wave can be accounted 
for by applying fhe theory of Reference 6. 

3.3.2 High Frequency Fluctuations: Acoustic Response - At plume fluctuation 
frequencies f > f^ , response time of the separated flow is slower than fluctuation of the 
plume. The separated region therefore does not respond quasi-statical ly to the plume motion. 
In the high frequency limit, f » f^ , the separation region may be considered to be 
steady (except, of course, for the pre-existing fluctuations which have been shown to be 
essentially independent of the plume), with plume fluctuations transmitted as acoustic 
disturbances. The fluctuating pressure in the separated region is then comprised of: 

• Fluctuation as described in Section 3.1 , modified according 
to Section 3.2. 

• Low frequency fluctuations where the separated flow responds 
quasl-statically to plume fluctuations, as discussed in Section 3.3.1 . 

• High frequency acoustic fluctuations superposed on the environment 
obtained from the first two above. 
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This last environment represents the jet noise of the rocket exhaust, and should be treated as 
a separate acoustic environment, rather than a part of plume-induced separation environments. 
The connection between this and separation is that the separation region provides a medium 
for upstream noise propagation at supersonic speeds. The model discussed so far does indicate 
that the acoustic environment may be divided into two regimes, which are worth pointing out 
at this time: 

• ^ or ^ se p^ f ^ fp| / fhe plume interface within the separated region may be 
considered to move solenoidally. The acoustic disturbance, within the 
context of the Lighthill theory (Reference 37), is then a distribution of 
synchronous sources along the plume interface. 

• For f > f | , the inhomogeneity of the plume must be accounted for. The 
sound source is then a distribution of quadrupoles. In this regime, the 
acoustic field must be treated as a full near-field jet noise problem. 
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4.0 APPLICATION TO SHUTTLE PRR CONFIGURATION 

4.1 Shuffle Flow Field 

Figure 28 shows the PRR shuttle configuration. It is comprised of four separate units linked 
together: an orbiter, an external tank, and two solid rocket boosters (SRB). The flow field 
for the mated launch vehicle will consist of the individual flow fields of each major body, 
with modifications due to interference between bodies. Because of the complexity of this 
vehicle, extensive theoretical and experimental studies will be required to accurately define 
the actual flow over the mated bodies. 

Despite the complexity of the interference regions of the shuttle flow field, there are some 
major areas where the flow will be simple enough so as to permit application of the simple 
axisymmetric plume model developed in this report. For zero angle of attack, the flow over 
the outboard sections of the SRB should be adequately described by locally axisymmetric flow 
about the centerline of each SRB. Provided the radii of the individual plumes are small enough 
so that they do not mutually interfere to a significant degree, the axisymmetric plume separation 
model developed in Section 2.0 may then be applied to the outboard section of the SRB with 
the same degree of confidence as the application to Saturn V separation presented earlier. 

To a lesser degree, the axisymmetric separation model may also be applied to the top surface 
of the orbiter. There is no direct interference in this region. However, the orbiter alone is 
not axisymmetric, the flow up to the flat top surface is over a relatively complex three-dimen- 
sional body, and the aft region is cluttered by the presence of the OMS/RCS pods and a 
vertical control surface. Therefore, although this calculation is performed below, the results 
should be considered to be an indication of trends rather than a quantitative prediction of 
separation point on the orbiter. 

4.1 .1 Plume -Induced Separation on SRB - To compute separation on the outboard 
side of the SRB, each SRB is treated as a single axisymmetric vehicle at zero angle of attack. 

The pertinent vehicle parameters for the SRB of the PRR configuration are summarized in Table 
I, based on data presented in References 38 and 39. 
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TABLE I 


SRB PARAMETERS 

Vehicle Radius 

81 in. - 2.06 m 

Thrust 

4, 162, 204 lb = 18.5 x 1 0 6 N 

y ex 

1.145 

Nozzle Area Ratio 

10 

Exhaust Mach Number 

3.18 

Chamber Pressure 

680 psia = 4.68 x 10 6 N/m 2 


Plume separation and fluctuation parameters were calculated for the due East trajectory used 
in Reference 7, shown in Figure 29. The results are shown in Figures 30 through 33. 

Figure 30 shows plume-induced separation length. Separation begins at about 95 seconds 
(altitude = 18 km f = 1 .8, q^ = 18,500 N/m 2 ), and encompasses the entire SRB at SRB 
cut-off of 145 seconds (altitude = 45 km, M^ = 4.6). Note that separation does not occur 
until well after conditions of maximum dynamic pressure, 82 seconds, altitude = 13 km, 

M co = 1 * 332 ' q co = 22,460 N//m2 ‘ C ° m Paring dynamic pressure, it may be concluded that, 
unless fluctuating pressure coefficients in the plume separation region are greater than 1 .33 
times those of flare-induced separation, plume-induced separation will not introduce more 

severe fluctuating pressure environments than those at the critical condition of maximum 
dynamic pressure. 


Figure 31 shows di f /(dp /q^ ). Using the feedback pressure coefficient shown in Figure 23, 
the change in separation length was obtained. Figure 32 shows < Lp > i , the increase in 
separation length due to the flexible plume boundary, along with boundary layer thickness at 
the separation point and root mean square displacement for an equivalent rigid flare. At the 
onset of separation, < Ai $ 2 > l is slightly less than < x| w > l, so that for perfect coherence 
shock oscillation distance is almost double that for a flare. Doubling oscillation distance 
would no more than double fluctuating pressure, or an increase of 6 dB. Adding 6 dB to the 
prediction presented in Reference 7 (and using the correct dynamic pressure for 95 seconds), 
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shock osci llaHon levels at the onset of separation will be approximately 161 dB. This is 3 dB 

higher than shock oscillation level of 158 dB predicted at maximum q in Reference 7. and 

00 

4 dB less than protuberance-induced homogeneous fluctuating pressure level of 165 dB at 

maximum q . 

oo 

Fluctuating pressure levels in the homogeneous separated region will be the same as for a rigid 
flare, hence will not be as great as those experienced at the critical condition of maximum 

q • 

00 

Figure 33 shows the change of separation length with fluctuation in engine chamber pressure. 
This prediction may be applied when data on SRB chamber pressure fluctuations are available. 

4.1 .2 Plume-Induced Separation on Orbiter - For the purpose of estimating 
separation trends, the present prediction scheme has been applied to the shuttle orbiter. The 
calculation was performed for an axisymmetric vehicle of radius comparable to orbiter 
dimensions. The vehicle and engine parameters used are listed in Table II. 


TABLE II 


ORBITER PARAMETERS 

Vehicle Radius 

1 22 in . = 3,11m 

Thrust 

1 .4 x 1 0 6 lb = 6.2 x 10 6 N 

r ex 

1.23 

Nozzle Area Ratio 

35 

Exhaust Mach Number 

4.05 

Chamber Pressure 

3000 psia = 2.065 x 10 ; N/m 2 


Figure 34 shows predicted separation length. Note that at 145 seconds, corresponding to SRB 
cut-off, the predicted separation point has not yet extended beyond the OMS/RCS pods and 
vertical control surface. This cluttered region is expected to experience p'otuberance-induced 
fluctuations and interference effects; therefore, plume induced separation will not make the 
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flow environment more severe. At later times, dynamic pressure is low enough that separated 
flow fluctuating pressures are not important. Therefore, it may be concluded that plume-induced 
separation due to the orbiter engines is not an important factor evaluating fluctuating flow 
environments. 


5.0 CONCLUSIONS 


The separated flow environment induced by underexpanded rocket plumes during boost phase 
of rocket vehicles has been investigated. A simple model for predicting the extent of 
separation was developed. The model employs a semi-empirical geometric representation of 
the plume interface. It was concluded that this simple representation, rather than a numerical 
method of characteristics computation as employed by other investigators, is sufficiently refined 
so as to be consistent with other approximations generally employed in the calculation of 
plume-induced separation. This simple model has the advantage of much greater economy of 
computation. Good agreement was found between present prediction of separation length 
and limited data from a Saturn V flight. 

The unsteady pressure field in plume-induced separated regions was investigated. Emphasis 
was placed on determining differences between fluctuations associated with plume-induced 
separation and those associated with an equivalent rigid flare. The frequency domain of 
possible fluctuating pressure/plume boundary coupling was investigated. It was found that 
coupling is limited to low frequencies, below the characteristic frequencies associated with 
the homogeneous separated flow fluctuations. It was further found that the main effect of the 
plume is to allow greater shock wave excursion distances than for an equivalent rigid flare. 
Differences in pressure fluctuations are limited to those associated with shock wave oscillation, 
for which the plume environment is more severe. The quantitative difference was found to be 
a function of vehicle and engine characteristics and flight conditions. The trend is for 
differences between the two cases to become less as separation length increases. 

The prediction scheme developed here has been applied to the PRR space shuttle configuration. 
The following results were found: 

• Separation due to the SRB engine begins at approximately 18 km, and 
extends over most of the vehicle by SRB cut-off (45 km). 

• Shock oscillation fluctuations associated with SRB plume-induced separation 
at the onset of separation are slightly less than 6 dB greater than those for an 
equivalent flare. At 18 km, the fluctuating pressure level is 161 dB; this 
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compares to maximum dynamic pressure (13 km) conditions of 158 dB 
for flare-induced shock oscillation and 165 dB for protuberance-induced 
separated flow. 


Shock oscillation levels exceed those for a rigid flare by lesser amounts 
as altitude increases above 18 km. 

Homogeneous separated flow pressure fluctuations induced by the SRB 
plume are comparable to those for a rigid flare. 

Separation due to the orbiter engines occurs at altitudes sufficiently high 
that it may be discounted as an important fluctuating pressure environment. 
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FIGURES 
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Figure 1 . Plume Induced Separated Flow 
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Figure 2. Schematic Representation of Shock-Wave/Laminar- 
Boundary-Layer Interaction, from Reference 11 
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Figure 6. Correlation of Plateau Pressure Ahead of a Compression Corner, from Reference 15 
as Collected in Reference 15, With Solid Points Due to Elfstrom 









Outer Plume Shock 



Figure 8. Inviscid High Speed Underexpanded Rocket Plume 
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Figure 9. Universal Plume Shape, from Reference 25 
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Figure 17. Typical Power Spectra in Region Ahead of a Supersonic 
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Three-Dimensional Protuberance, Showing Multiple 
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Figure 25. Separation Length Fluctuation, Boundary Layer Thickness at Separation, and Shock 
Oscillation Distance for Equivalent Rigid Flare, for Saturn V 
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Figure 27. Change of Separation Length with Change in Chamber Pressure 
for Saturn V 





28. Scnemaflc of PRR Parallel Burn Space Shuttle Launch Configuration 











Figure 30. Prediction of Plume-Induced Separation Point on Shuttle SRB 
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APPENDIX A 

DESCRIPTION Of 
COMPUTER PROGRAM 7J/003P-1 


WYLE LABORATORIES 
COMPUTER PROGRAM DESCRIPTION 


Program Number: 
Author: 

Date: 

Source Language: 
Computer: 


73/003P-1 
K. J. Plotkin 
February 1973 
FORTRAN IV-H 
XDS Sigma V 


1 .0 PROGRAM TITLE 
Plume Separation 

2.0 PURPOSE 

Given axisymmetric launch vehicle and trajectory parameters, the program computes 
properties relative to plume-induced separation. This includes separation length and 
plateau pressure, plus various properties of the plume and parameters describing the 
coupling to the unsteady pressure field. 

3.0 METHOD 


The analytic models incorporated in this program are described in detail in Wyle 
Labora 1 ori es Research Staff Report WR 73-3. 

4.0 COMPUTER CONFIGURATION 

The required hardware is: XDS Sigma V computer with 16 k core, card reader, and 
line printer. 

5.0 PUNCHED CARD INPUT 

Vehicle and trajectory parameters are input ?n metric units. Specific format for each 
card fs as follows: 


Card 1 - Descriptive title of vehicle. Up to 80 alphameric characters 


Card 2 


Columns 1-10 
Columns 1 1 -20 
Columns 21-30 
Columns 31 -40 
Columns 41 -50 
Columns 51 -60 


Vehicle radius, meters 

Total thrust, Newtons 

Plume drag to thrust ratio 

Exhaust gamma 

Exhaust Mach number 

Engine chamber pressure, N/m r; 


Format for each 


item on this card is 


FI 0.0. 
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Card 3 


Descriptive title of trajectory. Up to 80 alphameric characters 

N, the number of trajectory points to be read and computations 
performed for. Format 15. 

Cards 5 through 4+N — Trajectory parameters 

Columns 1-15: Time After Lift-Off, seconds 
Columns 16-30: Altitude, kilometers 
Columns 31-45: Flight Mach number 
Columns 46-60: Free-stream dynamic pressure, N/m 2 

Format for each item on these cards is E 15.10. 


6.0 OUTPUT 

The first items output are a table of atmospheric properties and a table of separation 
conditions contained in the program. This is followed by output of the vehicle 
parameters contained in Card 2, prefaced by the title input in Card 1 . Th it'e 
input in Card 3 is then printed, followed by a listing of the trajectory pare.neters 
input. Five (5) tables of computed properties are then printed. Each of these 
tables consists of the time after lift-off, followed by various properties. The 

Pr,n ^Li^ CienHy annofated 50 as be self-explanatory; the column headings 
are FORTRAN representations of various expressions defined in WR 73-3. A sample 
output list is presented herein. Written below each column is the exact notation 
used in WR 73-3, or a definition for those items not specifically denoted in the report. 

It should be noted that the program does not explicitly indicate whether or not the flow 
is separated. The output listing must be examined to determine this. The computed 
values represent separated flow when the plume angle is greater than the separation 

angle. In the sample output, flow is not separated at 84 seconds, but is at all other 
computation times. 
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